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We present an exact solution of the averaged Einstein's field equations in the presence of two 
real scalar fields and a component of dust with spherical symmetry. We suggest that the space-time 
found provides the characteristics required by a galactic model that could explain the supermassive 
central object and the dark matter halo at once, since one of the fields constitutes a central oscillaton 
surrounded by the dust and the other scalar field distributes far from the coordinate center and can 
be interpreted as a halo. We show the behavior of the rotation curves all along the background. 
Thus, the solution could be a first approximation of a "long exposition photograph" of a galaxy. 
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Doubtless, one of the most interesting open questions in physics by now is the one concerning the nature of the 
dark matter in the Universe. In a series of recent works we have proposed that the dark matter in Cosmos is of scalar 
field nature jl]-Q ■ Following an analogous procedure as in particle physics, we wrote a Lagrangian with all the terms 
needed to reproduce the observed Universe. In particular, using the scalar field potential for the scalar field dark 
matter (SFDM) 
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U($) = U [cosh(Ay^$)-l] (1) 



where k — 871- /M 2 ; , we were able to reproduce all the successes of the standard Lambda cold dark matter model 
( ACDM) at cosmological level very well including galactic scales [|I]j2] . The free parameters of the scalar potential 
Vo and A can be fitted by cosmological observations obtaining A ~ 20.23 and Vb — (3 x 10~ 27 Mpi) 4 , being Mpi the 
Q H . Planck mass ||. The mass of the scalar particle is then m$ ~ 9.1 x 10~ 52 Ajfp; = 1.1 x 10~ 23 eV (compare this value of 
the scalar mass with that in ||) . Under galactic scales there are some differences between the ACDM and the SFDM. 
bJQ' The self interaction of the scalar field $ could explain the observed dearth of dwarf galaxies and the smoothness of the 
galaxy core halos as well . The question whether this scalar field is from fundamental origin or it is a combination of 
• . some other fundamental particles, is open. Nevertheless, in order to reproduce the high resolution N-body numerical 
simulations of self interacting dark matter Q, we found that the renormalization scale of this scalar field is of the 
same order of the Planck mass 0. This would suggest that this scalar field dark matter could have a fundamental 
origin. We have also studied the scalar field hypothesis at galactic level in . The idea follows the standard idea of 
galaxy formation, namely that scalar field (dark matter) fluctuations are the responsible for the origin of the galaxies. 
In the case of the scalar field potential (1), we have cosh(A$) — > cosh(A(<i> + (5$)) ~ exp(A<5<i>) for regions where the 
scalar field fluctuation dominates <5$ > Of course, as in the standard theory of galaxy formation, the dark matter 
fluctuations are of different size in different regions of the Universe, for different galaxies. Therefore, at galactic level, 
we have a scalar field potential which depends on local variables. Thus, the exponential potential approximates the 
cosh potential in some regimes of the scalar field and we could develop some interesting aspects of a scalar dark matter 
halo in galaxies. 

On the other hand, new observations show that supermassive central objects lying in active galactic nuclei seem 
to be correlated with the velocity dispersion of the dark matter composing the dark halo, suggesting that the central 
object was formed at the same time than the halo ||. This is probably contrary to the standard idea about galactic 
nuclei that proposes the existence of a central black hole, since it should be formed a time after the disc, i.e. much 
more after the formation of the halo. Going further, it has been also shown that the supermassive objects which are 
supposed to be in the center of galactic nuclei, could be boson stars, obtaining the same predictions for the rate of 
accretion of matter due to the presence of a completely regular space-time background without surface or horizon || . 
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The question we are facing now is whether there is a dynamical mechanism that could provide a realistic scenario 
of galaxy formation using the scalar field dark matter hypothesis. First of all, a complete evolution of galactic and 
under galactic fluctuations belong to the non-linear regime of perturbations. The right answer would be provided by 
numerical evolution of Einstein's equations. Fortunately, a partial answer is given in numerical research on Einstein's 
equations developed since 1990. In particular, the collapse of a scalar field has been studied deeply in [jlO] ^2|, and 
it was found that there are final equilibrium and stable configurations for collapsed scalar field particles: boson stars 
(when the scalar field is complex) and oscillatons (when the scalar field is real and time-dependent), both of them 
being formed through a process called gravitational cooling Based on these ideas, we present in section [0] the 
motivations inviting us to build a galactic model with scalar fields, followed by the solution to the model in section 
III; physical features of the model are shown in section [ivl and finally we draw some conclusions. 



II. THE GALACTIC MODEL WITH SCALAR FIELD DARK MATTER 

We first recall the main results obtained in . Through the gravitational cooling process, a cosmological fluc- 

tuation of scalar field would collapse to form a compact oscillaton by ejecting part of the scalar field. This ejected part 
of the cooling process would carry the excess of kinetic energy out and can play the role of the halo of the collapsed 
object. The final configuration then consists of a central compact object, an oscillaton, surrounded by a diffuse cloud 
of scalar field, both formed at the same time due to the same collapse process. This would provide a correlation 
between the central object and the scalar halo. 

If the central object is an oscillaton, its formation is due to coherent scalar oscillations around the minimum of the 
scalar potential (Q). The scalar field <& and the metric coefficients (considering the spherically symmetric case) are 
time dependent and it was proved that this configuration is stable, non-singular and asymptotically flat ]l0[ ]. But, 
this could not be the final answer because, from the quantum-mechanical point of view, a Bose-Einstcin condensate 
should form and we must take into account that the scalar field <& can also decay or self annihilate [ p^p4| ] . In such 
a case, we should have to consider the scattering cross section <r$^$ in order to calculate the relaxation time of 
such condensation. In the models treated in |l3|-|l5[], the coefficient of <i> 4 is the responsible both for the size of the 
compact object and for the scattering cross section er$^$. With potential ([!]) it is necessary to take into account all 
the couplings 0. Nevertheless, it is possible to find the relaxation time smaller than the age of the Universe if [O 

g > l(T 1B (m»/eV) 7/2 (2) 

with g being the effective $ 4 -theory coupling. Then, considering the value m$ ~ 10~ 23 eV (see ||), g > 10~ 96 . This 
condition is well satisfied for our case since g ~ 10~ 48 0. This ensures that the relaxation time is shorter than the 
age of the Universe. 

In the case studied in Jll]], a massive scalar field without self- interaction collapses ejecting out 13% of the initial 
configuration of scalar field. The critical mass of the final configuration for the oscillaton depends on the mass of the 
boson JXOj ] . For our model (1), the mass of the scalar field is m — 1.1 x 10~ 23 eV and it is found that 

M 2 , 

M crit ~ 0.6—^ ~ 10 12 M Q (3) 
m 

which is a surprising result: the critical mass needed for the scalar field to collapse is of the same order of magnitude 
than the dark matter contents of a standard galaxy. Then, we expect that the fluctuations of this scalar field due to 
Jeans instabilities will in general collapse to form objects with mass of the same order than the mass of the halo of a 
typical galaxy. This is the reason by means of which we expect that the SFDM model can also work at galactic level. 
Summarizing, we consider two working hypothesis up to now. First, we identify the formation of a central compact 
object and a halo with the gravitational collapse of a scalar field. The compact object could be an oscillaton (since 
we are dealing only with a real scalar field) or a Bose-Einstein condensate. Second, we identify the ejected scalar field 
with the halo of this galaxy. 

In a realistic model the metric and fields should depend on time and a complete study would involve numerical 
calculations within and beyond General Relativity. One alternative is to study the behavior of the galaxy numerically 
with all the hypotheses stated above. This procedure has the advantage that we do not need to eliminate any a 
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priori consideration, but we can loose some important physical information inside of the numbers we obtain. Other 
alternative is to perform some approximations on the metric and fields and find exact solutions. These approximations 
could also lead us to loose crucial information of the system but we can keep physics under control. In our opinion, 
both procedures must be developed in order to have a better understanding of the hypothesis. 

In this work we will adopt the second alternative, i.e. we will build a toy model for case stated above in purely 
geometrical terms and considering only the final stage of the collapse. We let for a future work the dynamical evo- 
lution |T3]. To start with, wc support our toy model on the numerical results studied in [jlO] ^2|. First, since the 
time-dependence of the metric in an oscillaton is quite small |l0| , we suppose then that the center of this toy galaxy is 
an oscillaton which oscillates coherently but considering a static metric. This is an approximation because neither the 
galactic nuclei nor the oscillaton are expected to be static. However, for the purposes of this analytic work, we suppose 
that the dynamics of the oscillaton can be frozen in time in a way we explain below. Second, we do not expect the 
scalar halo to possess the same properties than the collapsed oscillaton; in some sense, they must be different. Thus, 
we will consider the scalar halo as another scalar field. Third, baryonic matter is considered to lie in the galaxy. We 
suppose that only the baryonic matter at the galaxy center and the bulge will contribute essentially to the curvature 
of the space-time of the galaxy. This baryonic component is assumed to be dust. As in previous works HQ, we let 
the luminous matter around the galaxy as test particles, i.e. they do not essentially contribute to the curvature of 
the space-time. 



III. THE ANALYTICAL SOLUTION 



The Lagrangian density we have to start with reads C = R + Cs fo + £s fh + £dust , being R the scalar curvature 
and the former terms correspond to the Lagrangian of the scalar field oscillaton (SFO), the scalar field halo (SFH) 
and dust, which should play the role of a galactic bulge. For the scalar objects, the following expressions are available 

Csfo = -\{d^d^) - V($), Csfh = -\{d^i>) - Uty) (4) 

Here $ is the scalar field corresponding to the oscillaton and V(4>) its potential of self-interaction, ip is the scalar 
field describing the scalar field halo and U(ip) its respective scalar potential. As a consequence of the non-coupling of 
the fields in the total Lagrangian, the Einstein's equations are written G M „ = no\T^° + T^ H + T*f s *], with 

= d^dv* - l/2 9fw [cr$d CT $ + 2V(<S>)} 



T^ H = d^d u il> - l/2 9llv [&>^d a <4) + 2t7(V0] 

-idust 



T^ st = du^u u 

being d the density of the dust and u a the four velocity of its particles. Since the scalar field $ is time-dependent 
oscillating coherently around the minimum of its scalar potential, we can write $ = P(r) cos(ivt) (see fl2| ). In 
order to handle Einstein's equations we average them during the period of a scalar oscillation, that is, we take 
< G>„ = kq\T^° + T^f H +T^ st ] >. This procedure gives the lowest order approximation for an oscillaton and we 
are left with time-independent differential equations [fL2| ; in the notation of reference p2| where the metric functions 
are written as g(r, t) = <?o( r ) + 9i( r ) cos(wt) + ... we are looking only for go terms of the metric and fields, which are 
the dominant ones. Besides of Einstein's equations there are two Klein-Gordon (KG) equations for the scalar fields, 



$; % + -7X = 0, + — =0 ( 5 ) 



The KG on the left depends on time, but it can be shown that one can factorize an overall cosine-term and then the 



resulting differential equation is time- independent 12 



Then, starting from a spherically symmetric space-time, using the harmonic maps ansatz we were able to find a 
solution of the system. In few words, the main idea behind the harmonic maps ansatz is the reparametrization of the 
metric functions with convenient auxiliary functions which will obey a generalization of the Laplace equation, along 
with some consistency relationships; the latter are usually quite difficult to fulfill, and great care and intuition must 
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be taken in order to get a system of equations both workable with and interesting enough |17[ . The exact solution of 
the averaged Einstein's field equations is 




(1 - v 2 a )KQ 




UW = ~ t^exp _,/ (6) 



for the scalar field ip. The scalar field <E> for the oscillaton we obtained is 

u 2f r 2 + h 2\-< r p2 /j + 2) r 2 + fe 2 A _ 2\ 

VY$M) = -1/4 i h 1/2 -i — — 

^ ^ 7 (-r + 2M)B (r 2 +& 2 ) Ko r 2 (l-v 2 ) 

- (v 2 + l) 2 r 4 + M«„ 2 (1 + 2 v 2 n ) r 3 -2b 2 (l + 2 v 2 ) r 2 + 5 « 2 Af6 2 r - 6 4 

| i/o \ Q / a V tty v ct/ a /y\ 

K r 2 (l-w4)( r 2 +6 2) 2 

where the function P is given up to quadratures 




P= I I 2 {l ~ Vt) 2 (1 ~ <] ^ - VIM (3 - 2 ^) r3 + 2 r2fe2 - 3 v * Mb2r + 54 ir f8) 
' u (r-2M)K r ( r 2 + b 2 ) 2 (r - 2 M) n a r 



being v a an asymptotic value of the tangential velocity of the test particles in our toy galaxy. The density of the dust 
is given by 

1 uj 2 (r 2 + b 2 y" l rP 2 

d= o-l/2 — - — -, ; 

K„r 2 ' B (r - 2 M) 

(l - vi) r A - v 2 a M (3-2 v 2 ) r 3 + 2 b 2 (l - v 2 a ) r 2 + v 2 a Mb 2 r + b 4 

K r 2 {l-vi) (r 2 +b 2 ) 2 (9) 

Finally, the corresponding line element reads 

ds 2 - -B (r 2 + b 2 Y« (l - 2 -^j dt 2 + {1 A0 m) dr2 + r 2 dtf (10) 

with dfl 2 = dO 2 + sin 2 8dp 2 and M is a constant with the interpretation discussed below. This metric is singular 
at r = 0, but it has an event horizon at r = 2M. This metric does not represent a black hole because it is not 
asymptotically flat. Nevertheless, for regions where r <C b but r > 2M the metric behaves like a Schwarzschild black 
hole. Inside of the horizon the pressure of the perfect fluid is not zero anymore, thus our toy model is valid only in 
regions outside of the horizon, where it could be an approximation of the galaxy. Metric ( |l0| ) is not asymptotically 
fiat, but it has a natural cut off when the dark matter density equals the intergalactic density as mentioned in ||. 



IV. PHYSICAL FEATURES OF THE MODEL 



The dust density and the oscillaton scalar field potential depend on the value of the function P. We find that 
P goes very rapidly to a positive constant value, which depends on the boundary conditions one imposes. Let its 
asymptotic value be Pq when r 3> 2M whose interpretation would be the asymptotic amplitude of tp. 

In order to understand the other parameters of the metric, let us proceed in the following way. It is believed that 
in a standard galaxy, the central object has a mass of M ~ 2 — 3 x 10 6 M Q ~some a.u. Far away from the center of 
the galaxy, say from Ipc up, the term 2M/r <C 1. In this limit metric ( |ic| ) becomes 

ds 2 = -B (r 2 + b 2 ) v «dt 2 + A Q dr 2 + r 2 dVt 2 . (11) 
This space-time is very similar to metric (18) of reference W, but now with the potential 
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U{ip{r)) = f° 2 j- (12) 

being both solutions the same in the limit r — ► oo. This implies that = 1 — 1>*, recovering in this way the asymptotic 
results shown in J|. Parameter 6 is related to parameter 6 of metric (21) in reference j§, where it acts as a gauge 
parameter. Of course this metric is only valid far away from the center of the galaxy. With parameter b it is now 
possible to fit quite well the rotation curves of spiral galaxies. Therefore metric (|l^) could not only represent the 
exterior part of the galaxy, but it could be a good approximation for the core part of it as well. Let us see this point. 

The rotation curves v rot seen by an observer at infinity for a spherically symmetric metric are given by v = 



y/ r 9tt,r/{^gu) @- For metric (10) such result reads 

v rot {r) 



v 2 a {r - 2M)r 2 + M(r 2 + b 2 ) 
(r - 2M)(r 2 + b 2 ) 



formula that allows one to fit observational curves. In Figure 1 it is shown the change of the rotation curve when 
the value of M changes, and it is evident that such fact would affect only the kinematics in the central parts of the 
galaxy, exactly in the same way as the mass of the central object should do M. 




FIG. 1. Rotation curves provided by the line element (Jlo|) , for different values of the mass of the central object: 
M = 3,5,7,9,11 x W 6 Mq. In the left side it is shown the rotation curve in the innermost part of the galactic nuclei, in 
the middle an intermediate radial distance an finally in the right hand side it is shown the properly called rotation curve 
associated to the dark matter. 

Nevertheless, it calls the attention in Figure 1 that there is a region where the velocity of test particles is near to 
zero, and in fact the rotation curve decays in a Keplerian way. Let us see what happens in such a region near the 
galactic center: the factor (r 2 + b 2 ) v ^ in g t t is almost 1 for r ~ lOpc for reasonable values b ~ kpc and v a ~ 10 -3 , 
therefore the time-like geodesies are determined by the factor 1 — 2M/r. Under such condition the contribution of 
-0 is accurately zero for an observer far away from the galaxy. The interesting situation comes after calculating the 
density of all the components to the central density, it reads 

, M(r-2M)(r 2 + b 2 ) v «~ 2 B , , , 

Pdust + Psfo + Psfh = -1/2 g 7- tt (-4 r v a 

r b K (1 - «*) 

+b 2 (4 - 7vl + 2 «*) r 4 + 3 6 4 (3 - v 2 a ) r 2 + 3 b 6 ) (13) 

where the approximation n a «l and (r 2 + b 2 ) v ^ ~ 1 is considered to be valid for reasonable values of the parameters. 
Moreover, we can infer that the total energy density as seen by an observer far away from the center of the galaxy 
diverges when r approaches to zero and goes rapidly to zero as 1/r 5 when r — > oo. But an observer at infinity measures 
the mass of the central object, it has not the possibility to distinguish between the density of each component. 

Now let us explore an ADM-like concept of mass associated to the central region; we stress that metric ( [uj| ) is not 
asymptotically flat and therefore the ADM concept of mass strictly fails to be valid, so we recall that there is a region 
into the interval 2M <C r < b geometrically almost flat, where we will define a useful infinity oof that should serve to 
calculate the mass of the central configuration through the standard metric 
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rlr 2 

ds 2 = -e 2S dt 2 + _ + (1 - a)r 2 dtf (14) 

V T ' 

where m = m(r) is interpreted as the mass function and 5 — S(r) as the gravitational potential. This form of the 
metric is convenient because in this coordinates fn r = 4irr 2 pT, where pt is the total density of the object. This 
interpretation is correct in regions where the space-time is almost flat, i.e. far away from the horizon r = 2M. 
Close to the horizon or inside of it, function m is a quantity that should be similar to the mass of the object, but 
it is not since it contains the contribution of all the components together; in this region, where the curvature of the 
space-time is huge, the volume element is different from 4jrr 2 dr. Furthermore, inside of the horizon we are not able 
to know the real physics of the object. On the other side, far away from the center of the toy galaxy, this function 
can be interpreted as the mass of an infinitesimal shell at radius r. Anyway we will call function m the mass function 
everywhere. Thus, the ADM-like mass is obtained by Madm = hm m. We perform the coordinate transformation 

r — >oo f 

\[A^r — ► r, \fA^b — > b in order to compare metrics (|o|) and (|l4|). We obtain 



ds 2 = ^ (r 2 + b 2 ) v ° 1 — - dt 2 + =— + — d9, 2 15 



with Aq — 1 — v^- Thus, Madm = \/l v iM. Probably an observer at any huge r would see the mass Madm at 
the center of the galaxy. 



There is a strong reason to consider with reserv the analysis near the horizon: we are considering a time averaged 
photograph of the space-time, which in the center would be strongly time dependent due to the presence of the 
oscillaton, and the stationarity fails to be a good approximation, thus it is not possible to know under the present 
approach which are the exact roles played by the binding energy of the oscillaton, the negative energy of the dust and 
the self-energy of the whole central system; such features could be known by performing the evolution of Einstein's 
equations and in fact represent a topic itself which will be discussed elsewhere |l6| . 

For this model, there is a contributor which does not fulfill the energy conditions, this is the reason to call it "toy 
model" , although an observer at infinity will see the sum of all energy densities of the components. In other words, 
the amount of matter, with negative or positive energy density, is the quantity which determines the value of m(r), 
the "balck hole mass" concentrated at radius r of the toy galaxy, not the contributors separately. Only if the observer 
could measure the contributions of the energy density very close to the center of the galaxy, it could recognize them. 
At infinity, the observer will only measure Madm, he. it will see a Black- Hole-like metric at the center of the galaxy 
which horizon lies at r = 2Madm- 

We present in Figure 2 the situation when parameter b takes different values, playing thus the role of the core radius 
in the usual dark halo hypothesis. Here we show that for an observer at infinity, the test particles close to the galaxy 
center behave as if there were a black hole of mass Madm m the center. The event horizon avoids that this observer 
could see inside of the horizon surface, the velocity of test particles is higher for closer test particles and the rotation 
velocity decays as 1/a/t in this region. For particles far away from the center, the rotation curves behave just as the 
rotation curves given by the contribution of the dark matter halo in a typical galaxy. 




FIG. 2. We show the rotational curves when the parameter 6=1,2, 3, 4, 5kpc. It is evident that in such case the shape of 
the curve associated with the dark matter changes, but in the center remains unchanged. 
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In Figure 3 the fit of the curves is done using the observed rotation curves of some dwarf galaxies, whose dark 
matter contribution is extremely dominating and therefore are considered as the test of fire for a dark matter model 
in galaxies. In general, for disc galaxies, the fit of the rotation curves using this metric is analogous to that in reference 
H . It seems then that metric ( |l0| ) is a good approximation for some late stadium of the space-time of a spiral galaxy; 
it is a good approximation of a "long exposition photograph" of a galaxy. 
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FIG. 3. Rotation curves of three dwarf galaxies are shown. The units are as usual: kms 
the horizontal one. 



in the vertical axis and Kpc in 



V. DISCUSSION 

The configuration we found is as follows: in the innermost region of the toy galaxy, it lies a time averaged oscillaton 
and dust, providing the picture of an object very similar to a black hole for an observer far away. Moreover, the 
resting scalar field should accommodate in the outer parts. This picture is precisely very similar to the one required 
for a current galaxy, which center would be dominated by a supermassive object and it should posses a halo region 
dominated by the presence of dark matter. The corresponding line element ( |Io| ) is free of naked singularities and is 
not asymptotically flat. Instead, it behaves exactly as a black hole near the coordinate center and it is asymptotically 
constructed to provide the flat curve condition Q . 

The main difference between this results and previous research supporting the scalar field dark matter hypothesis 
in galaxies, is precisely that in this case we have been able to extend the solution towards the galactic core. We have 
shown two important statements: First, an averaged spherically symmetric source field exists that supports a real 
scalar field under the assumption of non asymptotic flatness. Second, a fully relativistic treatment is able to explain 
the kinematic of test particles in the whole background space-time of a galaxy at once. 

Therefore, we believe that the metric presented here could be a first approximation of the galactic space-time 
provided the presence of coherent field (a Bose-Einstein condensate) instead of usual matter. Even when our solution 
accounts for negative energy densities of some components, observers at infinity are not able to know it. Furthermore, 
for such observers the test particles around the central object behave as if there would be a black hole of mass M^dm- 
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